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Abstract 

We have augmented the NASA DC-8 AIRSAR instrument with a pair of C-band antennas displaced across 
track to form an interferometer sensitive to topographic variations of the Earth’s surface. The antennas were 
developed by the Italian consortium Co.Ri.S.T. A., under contract to the Italian Space Agency (ASI), while 
the AIRSAR instrument and modifications to it supporting TOPSAR were sponsored by NASA. A new 
data processor was developed at JPL for producing the topographic maps, and a second processor was 
developed at Co.Ri.S.T.A. All the results presented below were processed at JPL. During the 1991 DC-8 
flight campaign, data were acquired over several sites in the United States and Europe, and topographic 
maps were produced from several of these flight lines. Analysis of the results indicate that statistical errors 
are in the 2-3 rn range for flat terrain and in the 4-5 m range for mountainous areas. 

Introduction 

We have implemented an interferometric synthetic aperture radar system for topographic mapping ap- 
plications (TOPSAR) on the NASA DC-8 aircraft. NASA /JPL currently operates a multi frequency (P, 
L, and C bands), multipolarimetric radar (AIRSAR) on board this aircraft. The TOPSAR implementa- 
tion uses much of the existing AIRSAR hardware, although several modifications were required to optimize 
performance in the topographic mapping mode. When in use, TOPSAR effectively replaces the C-band 
polar imeter instrument, but the remaining L- and P-band systems are undisturbed and operate together 
with the topographic mapper* therefore the combined instrument produces simultaneous L- and P-band 
fully polarimetric plus C-band VV polarization backscatter images in addition to the topographic product. 

In this paper we first describe our implementation, including system design parameters, required mod- 
ifications to the existing AIRSAR hardware, and new antenna design. We then discuss data processing 
strategy and display some example interferometric images. Finally, wo analyze an image acquired over Ft. 
Irwin, California, in terms of height accuracy. 

Background 

Interferometric radar has been been proposed and successfully demonstrated as a topographic mapping 
technique; by Graham [1], Xebker and Goldstein [2], and Gabriel and Goldstein [3]. Two distinct imple- 
mentation approaches have; been discussed for topographic radar interferometers which differ in how the 
interferometric baseline is formed. In one case a single; antenna and radar system illuminates a given surface 
at two different times but with nearly the same viewing geometry, forming a synthetic interferometer. This 
cast; lias been implemented both from spaceborne platforms (sec, for example, Goldstein et al. } [4]) and also 
from aircraft (Gray and Farris- Manning [5]). In the second approach the baseline is formed by two physical 
antennas which illuminate a given area on the ground simultaneously. This is the approach used originally 
by Graham [1] and also by Zcbker and Goldstein [2] lor the NASA CV-990 radar; it is the approach used 
here. 


Review of theory 

Consider a set of two antennas A1 and A2 as shown in figure 1. The surface topography is given by z(y), 
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h is the aircraft altitude, the baseline distance is B, the range to a point on the ground is p , the look angle 
0 , and the angle of the baseline with respect to horizontal is a. Radar echoes are transmitted from antenna 
A1 and received simultaneously at A1 and A2, thus the difference in path lengths is 6 , which depends on the 
baseline distance, baseline angle, range, look angle, and the height of the point z(y). The measured phase of 
the interferometer is directly proportional to this distance, with the constant of proportionality . A little 
algebra and geometry yield the following equations for height as a function of these parameters: 
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where <fi is the measured phase, and A is the wavelength. 

Differentiation of (1-3) yields the error in height estimate as a function of the error in phase estimate: 
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where a ^ and a $ are the standard deviations of height and phase, respectively. 

There exists an optimum baseline distance that minimizes the total height error [Rodriguez and Martin, 
6]; this optimum is quite broad and for high (> 10 dB) signal to noise ratio systems can be anywhere in the 
range of 0.2 - 0.8 of the critical value which may be calculated by ( Zebker and Villasenor [7]): 


B c = ^~ tan 0 

tty 


(5) 


where Ry is the slant range resolution. Note that this equation differs from that of Zebker and Villasenor by 
a factor of 2 cos 0 sin 0 as they were restricted to horizontal baselines, used the ground rather than slant range 
resolution, and assumed a single antenna in repeat-track configuration rather than two physical antennas as 
we use in TOPSAR. Another significant error source results from errors in knowledge of the aircraft attitude, 
the most important component of which is roll angle, which in our system translates directly to an error in 
look angle. Again differentiation yields 

a h = psm Oao (6) 

Thus, good motion compensation is a prerequisite of any practical system. 

TOPSAR design 

Several TOPSAR system parameters are summarized in Table 1. 


Table 1. TOPSAR radar system parameters 


Parameter TOPSAR value 


Peak power, watts 

1000 

Pulse rate, IIz 

600 nominal 

Pulse length, p sec 

5.0 

Antenna length, m 

1.6 

Antenna width, m 

0.11 

Antenna gain, dB 

25 

Range bandwidth, MIIz 

40 

Receiver noise temperature, K 

2100 

Antenna baseline, m 

2.58 

Baseline angle (a), deg 

62.77 


Given the baseline separation and the signal to noise ratio wc can then analyze performance of the 
interferometer. Figure 2 ( Li and Goldstein [8], Zebker and Villasenor [7]) gives the expected phase error as 
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a function of signal to noise ratio and number of looks for a radar interferometer. Thus we may expect a 
phase error of 3.3° if the data are processed to 16 looks assuming an SNR of 13 dB, which implies a statistical 
error of 1.50 m in height precision. Knowledge of the roll angle is also needed for accurate data reduction. 
At present our system can correct for attitude errors to the 0.015° level, which (6) implies a height error of 
1.9 m. 

Figure 3 shows a set of data acquired over Fort Irwin in the Mojave Desert in California. The data shown 
represent ground coordinate rectified topographic maps in which the color contour interval is 6 m, thus one 
trip around the color wheel corresponds to 96 m for our 16 entry color table. The spatial resolution is 10 m. 
The accuracy of the map obtained over Ft. Irwin is described in the next section. 

Verification procedure 

Verification of this mapping technique involves comparing the radar-derived height map with pre-existing 
high-resolution digital elevation models (DEMs) and characterizing the differences. Comparison was done 
on a data set from the Ft. Irwin area in California using the USGS 7.5’ DEM with 30 m grid spacing and 
7 m rms height error. Figures 4 shows a sample cut through the DEM and radar-generated terrain maps, 
and the difference between them is shown in the upper part of the figure. It is clear from examination of 
this figure that the error is dominated by differences in areas of high relief, and is caused both by errors in 
overall cross-track slope and by position errors. 
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1. TOPS A R geometry. Radar echoes are transmitted from an- 
tenna A1 and received simultaneously at A1 and A 2. The phase 
difference of the two echoes is proportional to the difference in 
path lengths 5 , which depends on the baseline distance B, base- 
line angle Ct, range p , look angle 6 , aircraft altitude h, and the 
height of the point z(y). 



TOPSAR SNR, dB 

2. Sensitivity of phase standard deviation to signal- to-noise ratio 
and number of looks in processor. Increasing number of looks is 
an effective means to reduce statistical variation, especially for 
the first 8 looks or so. 
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Along-track distance (m) 


3. Topographic data acquired over Fort Irwin in the Mojave 4. DEM / radar map height difference from the Fort Irwin data. 

Desert in California. The data shown represent ground coor- 
dinate rectified topographic maps in which the color contour 
interval is 6 m. 




